Numerical Simulation of 3D Acoustic Logging

A Finite-difference method for simulation of sonic waves propagating in a vicinity of a borehole filled with fluid and surrounded by a 3D heterogeneous elastic medium is developed. The method utilizes an explicit second order of approximation FD scheme on staggered grids that approximates the elastodynamic system of equations in cylindrical coordinates. A computational domain is surrounded by a special Perfectly Matched Layer for cylindrical coordinate system designed to attenuate waves reflected from outer boundaries. Parallelization is based on a domain decomposition approach and implemented with the help of MPI library. Results of numerical experiments are presented.